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To elucidate a novel pressure-temperature phase diagram of the quasi-one-dimensional mixed- 
stack charge-transfer (CT) complex TTF-CA, we study the quasi-one-dimensional spin-1 Blume- 
Emery-Griffith (BEG) model. In addition to the local charge transfer energy and the inter-stack 
polar (dipole-dipole) interaction, we take account of the inter-stack electrostriction effect. Using the 
self-consistent chain-mean-field theory, where the intra-stack degrees of freedom are exactly treated 
by the transfer-matrix method, we reproduce the gas-liquid-solid like phase diagram corresponding 
to the neutral (N), paraelectric ionic (I pM a), and ferroelectric ionic (If er ro) phases, respectively. We 
also give an explanation on the experimentally observed multicritical behavior and concomitant 
discontinuous inter-stack lattice contraction in TTF-CA. 



PACS numbers: 

I. INTRODUCTOIN 

In the "critical phase control technology," condensed 
molecular materials play quite a promising role, because 
molecular orbitals and stacking architecture are manip- 
ulable in a desirable way. To elucidate interrelation of 
constituent molecular structures and emergence of var- 
ious thermodynamic phases such as superconductivity, 
magnetism, and ferroelectricity is of great interest there. 
A neutral-to-ionic phase transition (NIT) in quasi-one- 
dimensional charge-transfer (CT) complexes comprising 
mixed-stack architecture of electron donor (D) and ac- 
cepter (A) molecules^ has played a key role in this field. 

In particular, phase control by pressure^ or laser 
radiation^ in the tetrathiafulvalene-p-chloranil (TTF- 
CA), that exhibits the NIT around 80K at ambient pres- 
sure, has attracted a great deal of interest. Very recently, 
Collet at al., 5 using highly refined time-resolved X- ray 
diffraction technique, have reported direct observation 
of a photo-induced paraelectric-to-ferroelectric structural 
order in the crystal. In the ionic phase, the D + A _ pair 
forms a dimer due to the electrostatic instability^ or sub- 
sequent spin-Pcicrls instability^ The ionized dimer on 
the DA chain carries a local electric dipole moment p 
with opposite directions depending on the dimerization 
patterns D + A ~ or A~D + . Once p acquires a macroscopic 
mean value r\ = (p) =/= 0, a spontaneous inversion sym- 
metry breaking (SISB) occurs and the system undergoes 
a phase transition to a ferroelectric-ionic (Iferro) phase. 
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The ionic phase itself is simply described by ionicity con- 
densation c = (p 2 ) ~ 1. Since r] is a symmetry-breaking 
order parameter but c is not, we expect that rj and c play 
separate roles. The appearance of two distinct order pa- 
rameters -q and c is a direct consequence of the degeneracy 
of the two configurations of dimerization pattern, I A (■ ■ • 
D+A - D+A - • • • ) and (• • • A-D+A_\D+ • • • )■ 

Recently, the respective roles of c and 77 have been high- 
lighted in both equilibrium^ and noncquilibrium^ pro- 
cesses. Using the neutron diffraction along with nuclear- 
quadrupole-resonance (NQR) measurements, Lemee- 
Cailleau et al. 2 * found a novel phase where the system 
is ionic but dipoles remain disordered, i.e., a paraelec- 
tric ionic (Ipara) phase. They proposed a pressure- 
temperature phase diagram of TTF-CA, where the N, 
Ipara, and If Grr o phases are like gas, liquid, and solid 
phases, respectively. The ferroelectric order is well sig- 
naled by the appearance of (0, 2fe+l, 0) Bragg peaks that 
indicate the inversion symmetry breaking. The "sublima- 
tion" line separating the N and If er ro phases continues 
up to a triple point (P t ,T t ) - (500MPa, 210K). Above 
the triple point, in addition to the "crystallization (or 
melting)" line, there appears a "condensation" line sep- 
arating the Ipara and If err o phases accompanied by a con- 
comitant discontinuous change of c, ending at a critical 
point (P C ,T C ) ~ (700MPa,250K). The purpose of this 
paper is to give a qualitative understanding of this phase 
diagram. This phase diagram, in addition to the phe- 
nomenological description offered in ref£. 

Since the SISB is prohibited by thermal fluctuations 
in a purely one-dimensional stack, inter-stack coupling 
is required to realize the SISB. In addition, the experi- 
mental observation strongly indicates that electronic and 
lattice degrees of freedom are coupled with each other in 
a unique manner. That is to say, upon crossing the tran- 
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sition lines in the gas-liquid-solid like phase diagram, 2,3 
the unit cell parameter b (for the axis perpendicular to 
the stack) exhibits about 0.5 % discontinuous contrac- 
tion at the condensation transition but exhibits continu- 
ous contraction at the crystallization transition. On the 
other hand, the unit cell parameter a (for the stacking 
axis) exhibits only continuous contraction at the conden- 
sation transition and below it remains almost constant. 3 
Now we are ready to ask the question: (1) what kind of 
inter-stack interactions are responsible for the occurrence 
of the Iferro phase, and (2) how the lattice anomalies are 
coupled to the phase transitions ? 

As for the first question, Luty et al.£& stressed that the 
inter-stack non-polar coupling 7 alone cannot drive the 
ferroelectric ordering and the dipolar coupling plays an 
essential role. As for the second question, Kawamoto et 
al& took account of the charge distribution on the atoms 
inside each molecule by an ab initio quantum chemi- 
cal method and elucidated the importance of inter-stack 
Coulomb attraction ~ — 0.14cV, which may cause inter- 
stack electrostriction (Coulomb-lattice coupling). 

II. QUASI-ONE-DIMENSIONAL 
BLUME-EMERY-GRIFFITH MODEL AND 
INTERCHAIN MEAN FIELD THEORY 

Now we shall set up a model. The ground-state energy 
of the mixed stacks has three minima as a function of the 
dimerization displacement, i.e., the N and the degenerate 
IA and IB states. The three states may be described by 
the spin-1 Ising variable pij = 0, ±1 on the ith dimer 
inside the jth stackiS* 3 ^ The charge transfer energy (A), 
the intra- (with subscript ||) and inter- (with subscript 
_L) stack dipolar (J) and non-polar (K) interactions, and 
the coupling with the electric field (E) are described by 
the quasi one-dimensoinal (Q1D) Blume-Emery-Griffith 
(BEG) models H = H\\+H±, where 

id 

-Ap*,. - Ejhj] , (1) 

h± = - [j± Pi , jPi ,j +1 + k ± p1 3 p1 3+1 ] . (2) 

The intra-stack dipolar interaction Jii is caused by 
coupling between the charge transfer and the lattice 
distortion^ while the inter-stack dipolar interaction is 
regarded as a direct interaction between the induced 
dipoles on adjacent stacks. The intra-stack couplings 
are much stronger than the inter-stack couplings, and 
the electric dipoles are aligned along the stacks. The 
energy cost to create one D + A~ pair is given in the 
limit of no molecular overlap by A = I — A — aV, 
where / and A denote the donor's ionization energy and 
the accepter's affinity, respectively, and aV denotes the 
Madelung energy^i Generally speaking, increasing pres- 
sure decreases the lattice spacing a and consequently in- 



creases V. Therefore, A decreases upon applying pres- 
sure. 

We treat the Hamiltonian by using the self- 
consistent chain mean-field theorjii^. Introducing the 
thermal averages, rj = (pij) and c = (pj ■), we have the 
effective ID BEG model, ' 

n f = -J2i J \\PM+i +K \\PiPi+i- K Pi-zpi] 

i 

+ Z -±NJ ±V 2 + ^NK ± c 2 , (3) 

where A = A — z±Kj_c and E = E — z±J±rj, with z± = 2 
being the inter-stack coordination number. Treating 7ijj ff 
exactly by the transfer matrix method, we obtain the free 
energy per site, 

/beg(?7,c,T) = -Tln\(A,E,T) + J ±V 2 + K ± c 2 , (4) 

where A(A, E, T) is the maximum eigenvalue of the trans- 
fer matrix for Wjj ff , given by 

/ e /3(J||+K||-A-B) e -/3(A+B)/2 e /3(- -A) \ 
T= e -/3(A+E)/2 j e -f}(A-E)/2 1(5) 

I e /3(-Jn+Kn-A) e -P{K-E)/2 e 0(J ]]+ K^K+E)J 

with j3 = 1/T. Possible phase diagrams of the BEG 
model have been extensively studied through mean- 
field theories^2*i£ renormalization-group^ and transfer- 
matrix methods^ For the parameter regions relevant to 
the present case, J||, Jj_, Ku, K±, and A are all positive, 
so that a solid-liquid-gas type phase diagram with proper 
slopes of transition lines is not obtained. 

III. INTERCHAIN ELECTROSTRICTION 

Then, we consider the inter-stack lattice degrees of 
freedom that have not explicitly been taken into account 
in (PJ. It is well known that an electrostriction effect 
potentially converts a continuous transition to a discon- 
tinuous one, since this gives rise to an additional neg- 
ative free-energy term that contains the forth power of 
the relevant order parameter^ In the present case, we 
phcncmcnologically introduce an additional free energy, 

c 2 1 

/eist(c,y) = -— K-fcy 2 , (6) 

oo + y 2 

where the first and second terms represent Coulomb at- 
traction between the nearest neighbor stacks^ and the 
elastic energy for the distortion in the inter-stack direc- 
tion. Note that / c ist(c, 0) has already been absorbed into 
K±. The lattice constant without distortion is bo, and 
y denotes the distortion. By minimizing / c i s t(c, y) with 
respect to y, we obtain the optimized lattice constant, 

b{T) = b + y(T)~b -2e clst b 2 c 2 1 (7) 
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where e c \ s t — 1/(2A;6q) is a small constant. We thus have 
the energy gain due to the lattice distortion, 



other, and the transition simply changes from discontin- 
uous to continuous at some critical value of AiiSii&ilii^ 



/eist(c, T) 



-EelstC 



(8) 



Now, solving the self-consistent equations is reduced to 
searching (c, rj) that gives the absolute minimum of the 
total free energy, f(rj,c,T) = /beg(»7, c i T ) + /cist(c, T). 
The ionic phase is characterized by the ionicity condensa- 
tion c = 1, while the ferroelectric phase is characterized 
byry^O. 

Note that, in the present scheme, any phase with 1 
is regarded as "neutral" and that the neutral phase has 
always r\ = 0. In the BEG model, because of three 
states pij = 0, ±1, c approaches the universal constant 
c = 2/3 in the high temperature limit, where the entropy 
term dominates the internal energy term. Therefore, 
in the parameter region where c continuously increases 
upon decreasing temperature, we have 2/3 < c < 1. 
In the experiments^ the ionicity continuously increases 
upon decreasing temperature and jumps from c ~ 0.3 
to c ~ 0.6 at the NIT. Thus, concerning the quantita- 
tive magnitude of the ionicity, there arises a difference 
between the experimental result and the present analy- 
sis. This apparent difference comes from the fact that we 
mapped the intra-stack CT transfer and the DA dimer- 
ization onto the simple spin-1 Ising variables. Therefore, 
we should regard the difference as an artifact of the clas- 
sical BEG model. 



IV. NUMERICAL RESULTS AND PHASE 
DIAGRAM 

From now on, we set J11 = 1 as an energy unit and 
the electric field E is set to be zero. In Fig. ^ we 
show the temperature dependence of c and r\ for var- 
ious magnitudes of A with K\\ — 0.4, K±_ — 0.06, 
J± = 0.03, and e e ist = 0.0095. We introduce the conden- 
sation temperature T con< i and the crystallization temper- 
ature T cryst . The ionicity jumps into c = 1 at T con d, while 
the ferroelectric order parameter acquires a finite mag- 
nitude i] / at T cryst . The ground state becomes ionic 
for A < 1.49. Both c and n exhibits a discontinuous 
change at the same transition temperature (the sublima- 
tion temperature) for 1.42 < A < 1.49. That is to say, 
Tcond = ^cryst. For A < 1.42, there appears a region, 
Tcryst < T < Tcond, where the system is ionic but still 
paraelectric. This region is identified with the / par a phase 
that is observed in TTF-CA under pressure. The point 
(At = 1.42, T t = 0.45) is identified with the triple point 
(indicated by "TP" in Fig.0). For A < 1.42, c still ex- 
hibits a discontinuous change at T con( j, but n continuosly 
evolves at T cryst , as shown in Fig. ^c). The discontinu- 
ity jump of c becomes smaller and seems to vanish, as A 
decreases, as shown in Fig.^d). We stress that this dis- 
continuity is a direct consequence of the weak but finite 
electrostriction effect. Without the electrostriction, as A 
decreases, T con d and T crys t continue to coincide with each 
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FIG. 1: Temperature dependence of c and n for various mag- 
nitudes of A with K\\ = 0.4, K± = 0.06, J± = 0.03, and 
e e ist = 0.0095. The condensation and crystallization temper- 
ature, T con d and T crys t, respectively, are indicated. Locations 
of the A values in (a)-(d) are indicated in the phase diagram 
of Fig. 

As clearly seen from {7|, about 0.5% discontinuous 
contraction of the inter-stack lattice constant (unit cell 
parameter b) is accompanied by the discontinuous jump 
of the ionicity. In Fig. [3 setting &o as a length unit, we 
show the temperature dependence of the unit cell param- 
eter b given by , using the same parameter set as that 
in Fig. n Although the magnitude of the discontinuous 
contraction depends on the parameter choice of £ e i s t and 
&o, the qualitative nature (b jumps at T com d) does not 
change. 

To see the discontinuity of the ionicity more closely, 
we show in Fig. the A dependence of the discontinuity 
at the condensation temperature, Ac. It is clearly seen 
that Ac decreases as A decreases and eventually reaches 
zero at A = 1.25. For A < 1.25, the condensation oc- 
curs without ionicity jump. Then, the lattice contraction 
at Tcond also becomes continuous. Therefore, A = 1.25 
with the corresponding T con d = 0.76 is idendentified with 
a critical point (indicated by "CP" in Fig.[3J). This re- 
sult is well consistent with the experimental fact that the 
ionicity jump finishes at the critical point. 2 

The dielectric constant is given by e = 1 + 
47ra, where the uniform polarizability is a = 

T^id [(PijPi,m> - (Pi,j)(pi,m)] = (c ~ ^/T. In 

Fig. 01 we show the temperature dependence of a for 
various magnitudes of A. It is seen that along the N- 
Ifcrro boundary, the polarizability exhibits a sharp single 
cusp at T con d = T cr y S t. For A c < A < A t , a discontinu- 
ous jump occcurs at T = T con d and a cusp at T = T C r ys t- 
The discontinuity at T = T con d finishes at A = A c . 
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In Fig. |SJ we show the phase diagram of the system 
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FIG. 2: Temperature dependence of the unit cell parameter 
b, b(T). Locations of the A values in (a)-(d) are indicated in 
the phase diagram of Fig. 
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0.0095. Regarding the decreasing A as increasing pres- 
sure, this phase diagram is consistent with the experi- 
mentally found, pressure-temperature phase diagram of 
TTF-CAiS The tiriple point, the critical point, and the 
observed inter-stack lattice contraction are reproduced. 
For simplicity, we here ignored the change of A due to 
thermal lattice contraction. Exactly speaking, to con- 
vert our A-T phase diagram to a P-T diagram, we need 
to take account of the temperature dependence of A, 
A(T). By appropriately treating A(T), we may ob- 
tain the corresponding P — T phase diagram satisfying 
the Clausius-Clapeyron relation. We stress that, even 
when we take this simple view, a qualitative nature of 
the phase diagram is not changed. Identifying the triple 
point (Af,T t ) = (1.42,0.45) with the experimentally ob- 
tained one (P t ,T t ) ~ (500MPa, 210K), we sec that our 
parameter choice here corresponds to K± = 28K and 
J ± = UK. 

Lajzerowicz and Sivardiere 1 ^ extensively developed a 
mean-field analysis of the BEG model and obtained 
liquid-gas-solid like phase diagrams on the P-T plane. 
However, they considered a lattice gas analogue of a sim- 
ple fluid, where the physical pressure of the lattice gas is 
simply given by — /, with / being the Helmholtz free en- 
ergy per volume. In the present context, the pressure of 
the spin system has no physical meaning and the phase 
diagram obtained by Lajzerowicz and Sivardiere cannot 
be applied to TTF-CA. 



FIG. 3: Temperature dependence of the discontinuity of the 
ionicity, Ac. 



CONCLUDING REMARKS 
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FIG. 4: Temperature dependence of the polarizability, a. 
Locations of the A values in (a)-(d) are indicated in the phase 
diagram of Fig. |K| 



In this paper, we showed that the inter-stack elec- 
trostriction causes the — c 4 term with a small coefficient 
and makes the phase diagram rich. We thus conclude 
that the inter-stack polar interaction together with the 
inter-stack electrostriction drives the discontinuous inter- 
stack lattice contraction and the multicritical behavior 
observed in TTF-CAi^ Our considerations are limited 
to classical models, i.e., we considered classical effective 
models, where all the microscopic (electronic or phonon) 
degrees of freedom are implicitly integrated out. To ob- 
tain the truly microscopic mechanism behind the mul- 
ticriticality, we need to get back to such a microscopic 
Hamiltonian as an extended Peierls-Hubbard models 
We would keep this issue for future study. 
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FIG. 5: Phase diagram for K\\ = 0.4, K± = 0.06, J ± = 0.03, 
and e e ist = 0.0095. The solid and dashed lines represent dis- 
continuous and continuous transitions, respectively. TP and 
CP represent the triple point and the critical point, respec- 
tively. Locations of the A values used in (a)-(d) of Figs. [5] 
and 0| are indicated by the horizontal arrows. 
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